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EQUILIBRIUM  NOISE  IN  ION  SELECTIVE  FIELD  EFFECT  TRANSISTORS 
Andre  Haemmerli,  Jtr{  Janata,  and  James  J.  Brophy 
Department  of  Bioengineering 
University  of  Utah 
Salt  Lake  City,  Utah  84112  U.S.A. 

ABSTRACT 

Noise  spectra  of  the  drain  current  of  ion  selective  field  effect 
transistors  (ISFETs)  in  equilibrium  were  measured.  The  spectra  differ  from  one 
type  of  ISFET  to  the  other.  The  ionic  activity  in  solution  does  not  modify 
perceptibly  the  noise  level  for  any  of  the  various  types  of  ISFETs.  Computer 
assisted  modeling  is  used  to  fit  experimental  data  and  to  assign  an  equivalent 
electrical  circuit  to  the  electrochemical  system.  The  overall  response  time  of 
the  system  is  shown  to  depend  on  the  input  capacitance  of  the  electrometer. 
INTRODUCTION 

A  considerable  amount  of  research  has  been  devoted  to  the  ion  selective 
field  effect  transistor  (ISFET), and  the  subject  has  been  recently  reviewed 
[ 1 ] .  Most  of  the  work  has  been  directed  to  the  description  of  the 
deterministic  response  to  activity  and  potential  changes  of  ISFETs  selective 
to  different  kinds  of  ions.  Fluctuation  studies  on  many  different 
electrochemical  systems  have  been  reported  in  the  literature.  Such  studies 
include  redox  reactions  at  metal  and  semiconductor  electrodes  [2-5],  the 
dissociation  process  of  beryllium  sulfate  [6],  the  nucl'  '  -»  and  crystal 

growth  of  metals  and  semiconductors  and  the  nucleation  ot  ).  j  in  lipid 
bilayers  as  well  as  the  kinetics  of  electrocatalysis  [7],  and  the  fluctuation 


analysis  in  biological  membranes  [8-10]  to  name  only  a  few.  These  studies  show 
that  the  measurement  of  random  electrical  fluctuations  (commonly  referred  to 
as  noise)  occurring  in  electrochemical  systems  can  lead  to  kinetic 
information,  some  of  which  is  not  accessible  by  deterministic  measurements. 

This  article  presents  a  study  of  random  fluctuations  in  the  drain 
current  of  ISFETs  in  the  range  of  1  Hz  to  25  kHz  for  several  different  ion 
selective  membranes.  ISFETs  are  normally  used  as  potentiometric  sensors,  which 
means  that  no  net  current  flows  through  the  membrane/solution  interface.  In 
most  ion  selective  devices  the  zero  net  current  condition  is  synonymous  with 
thermodynamic  equilibrium.  It  was,  therefore,  decided  first  to  investigate  the 
phenomenon  in  ISFETs  in  this  state.  The  overall  noise  exhibited  by  the  drain 
current  of  an  ISFET  is  the  sum  of  several  contributions.  These  contributions 
can  be  separated  as  follows: 

a.  Noise  generated  in  the  "solid  state"  portion  of  the  device  which 
includes  thermal  noise,  generation-recombination  noise  in  the  space  charge 
region  at  the  substrate-channel  interface,  and  1/f  noise  in  the  channel  [11, 
12]; 

b.  Electrochemical  noise  associated  with  the  ion/membrane  interactions 
which  can  exhibit  thermal  noise,  shot  noise,  and  1/f  noise; 

c.  Noise  generated  in  the  solution  and  at  the  reference  electrode,  as 
well  as  noise  resulting  from  the  fluctuations  of  the  biasing  elements  (power 
supplies,  resistances).  The  noise  spectra  for  several  kinds  of  ISFETs  were 
measured  and  interpreted  in  terms  of  the  electrochemical  properties  of  the 
membranes. 

NOISE  CHARACTERIZATION 

The  drain  current  1^(0  in  the  channel  of  ISFETs  can  be  expressed  by 


Equation  1: 


V‘>  =  :D  .  Kt) 
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Where  1^  is  the  average  value  and  i(t)  is  the  random  fluctuation  of  the  drain 
current.  The  random  fluctuation  or  noise  i(t)  can  be  quantified  and 
characterized  by  means  of  spectral  density  analysis,  correlation  analysis  and 
integral  spectrum  analysis  [10,  13,  14],  From  an  analysis  point  of  view,  there 
is  little  difference  between  methods, but  there  may  be  some  practical 
difference  which  favors  one  or  the  other  in  particular  circumstances. 
Furthermore,  information  obtained  with  one  method  can  be  transformed  into 
information  obtained  with  another  method. 

For  our  application,  we  used  spectral  density  analysis  mainly  because 
it  is  the  type  of  analysis  that  has  been  used  most  often  for  the  measurement 
of  noise  in  membranes  and  in  electronic  devices.  In  order  to  define  the 
spectral  density  of  i(t),  we  use  Rayleigh's  theorem  which  states  that  the  area 
under  the  square  modulus  of  a  function  is  equal  to  the  area  under  the  square 
modulus  of  its  Fourier  transform.  This  is  represented  by  Equation  2: 


r  +  °o  r  +  oo 

j  U<t)|2dt  =  j  |Kf)| 


df 
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Where  1(f)  is  the  Fourier  transform  of  i(t)  and  is  by  definition  given  by 
Equation  3: 
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where  w  =  2irf,  with  f  being  the  frequency.  Since  i(t)  is  a  real  signal  which 
exists  only  for  0  <  t  <  T,  Rayleigh's  theorem  can  be  rewritten  as: 


|  T  i2(t)  dt  =  |  | i C f )  |2  df 


(4) 
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The  average  value  of  i(t)  is  zero  and  thus,  the  average  value  of  i  (t)  is  the 
2  — 

variance  a  .  On  the  other  hand,  |I(f)|  is  an  even  function  of  frequency  which 
means  that  the  integral  over  all  frequencies  may  be  expressed  as  twice  the 
integral  over  the  positive  frequency  axis.  So,  using  Equation  4,  we  can  write: 
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jT  i2(t)  dt  =  |  |  |I(f)|2  df 
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By  definition,  the  spectral  density  S^f)  of  i(t)  is  given  by: 


(5) 


S  (f)  =  lira  |  |I(f)|  2  (6) 

1  x  00 

where  the  subscript  i  refers  to  the  spectral  density  of  the  current. 
Combining  Equation  5  and  6  leads  to: 
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This  shows  that  the  variance  of  a  signal  in  a  narrow  frequency  band  df  is 
given  by  S^fldf  which  clarifies  the  meaning  of  the  spectral  density. 

If  the  signal  is  deterministic  and  is  defined  by  an  analytic  expression, 
the  spectral  density  and  the  variance  are  absolutely  determined.  On  the  other 
hand,  if  the  signal  is  probabilistic,  such  as  random  noise,  the  spectral 
density  and  the  variance  will  be  different  for  each  time  interval  T.  However, 
if  the  random  signal  is  stationary,  which  means  that  the  variance  is 


independent  of  time,  the  different  estimates  for  spectral  density  and  variance 
will  tend  toward  a  mean  value. 


There  are  many  experimental  ways  to  estimate  the  spectral  density.  The 
common  feature  of  the  different  techniques  is  filtering  of  the  signal.  For  our 
measurements,  we  used  a  digital  spectrum  analyzer  (HP  3582A)  in  which  the  fast 
Fourier  transform  (FFT)  implements  a  set  of  parallel  filters  digitally  and 
basically  calculates  the  spectral  density  as  a  function  of  the  center 
frequency  of  the  filter.  For  more  information  about  the  characterization  and 
the  analysis  of  noise,  see  for  example,  DeFelice  [10]  or  Van  der  Ziel  [13, 

14]. 

EXPERIMENTAL 

Devices 

The  fabrication  and  modes  of  operation  of  ISFETs  have  been  described 
previously  [l].  Basically  an  ISFET  is  obtained  by  replacing  the  metal  gate  of 
the  insulated  gate  field  effect  transistor  (IGFET)  with  a  reference 
electrode/electrolyte/ion  selective  membrane  electrochemical  system. 

Figure  1  is  a  photograph  of  the  integrated  circuit  we  used  in  our 
experiments  and  which  has  been  described  in  Reference  1.  This  integrated 
circuit  is  composed  of  two  ISFETs  (Q^,  Q^)  and  two  IGFETs  (Q2»  Q^).  All 
devices  are  n-channel  depletion  mode  field  effect  transistors. 

The  devices  investigated  in  this  study  can  be  divided  into  four 
categories,  depending  on  their  gate  structure.  The  first  is  the  standard  IGFET 
structure,  which  is  considered  to  be  the  basic  device  against  which  ISFET 
performance  is  measured.  The  second  is  the  bare  gate  silicon  nitride  ISFET 
which  is  known  to  be  pH  sensitive  [15].  The  third  group  includes  ISFETs  with 


-5- 


I 


homogeneous  polymeric  membranes;  Ca  +,  K+,  and  Na+  selective  membranes  were 
used  [16-18].  We  have  also  measured  ISFETs  with  gold  covered  gate  which  are 
sensitive  to  redox  couples  (redox  ISFETs)  and  which  can  be  operated  either  as 
an  IGFET  or  as  a  redox  ISFET. 

Power  Spectrum  Measurements 

The  experimental  arrangement  used  for  electrical  noise  measurements  is 
shown  schematically  in  Figure  2.  The  gate  bias  was  provided  by  a  9  V  alkaline 
cell  (Ray-O-Vac  No.  A1604)  appropriately  filtered  by  the  low-pass  combination 
R^  and  C  in  order  to  avoid  possible  contamination  with  battery  noise.  The  -3dB 
frequency  f^p  for  this  low-pass  filter  is  given  by  Equation  8: 

f  -  - 1 _ 

LP  2  ir  R3  C  (8) 

which  with  our  values  gave  a  frequency  of  1.59  Hz.  For  measurements  on  ISFETs, 

the  reference  electrode  used  was  a  custom  made  Ag/AgCl  (saturated  KC1)  cell  of 

2 

surface  area  approximately  5  cm  .  The  drain  bias  was  such  that  the  devices 
were  operated  in  the  saturation  region.  The  voltage  (Figure  1)  consisted  of 
a  series  combination  of  four  1.5  V  alkaline  cells  (Duracell  MN  1300). 

The  drain  current  Ip,  which  was  maintained  at  100  yA  for  all  experiments 
developed  a  voltage  across  R^  (metal  film  resistor,  Dale  model  RN60D,  CMF60D) 
and  the  fluctuations  of  this  voltage  were  amplified  using  a  low  noise  pre¬ 
amplifier  (PARC  113).  The  gain  selected  for  all  experiments  was  1000  except 
for  amplifier  and  resistor  noise  measurements  where  it  was  10,000.  The  pre¬ 
amplifier  was  ac  coupled  and  its  -3dB  high  frequency  roll-off  and  low 
frequency  roll-off  were  300  kHz  and  0.03  Hz,  respectively.  The  amplified 
signal  was  analyzed  using  a  digital  spectrum  analyzer  (HP  3582A).  This 
instrument  has  a  frequency  range  of  0.02  Hz  to  25  kHz.  The  combination  of  the 
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pre-amplifier,  the  gate  bias  low-pass  filter  and  the  spectrum  analyzer  thus 
allowed  us  to  investigate  the  frequency  range  of  1  Hz  to  25  kHz. 

The  IGFET  (or  ISFET)  and  its  biasing  components  were  enclosed  in  a 
Faraday  cage  to  minimize  external  electrical  interference.  For  each  device, 
spectra  were  measured  in  the  frequency  range  1-10  Hz,  1-100  Hz,  1  Hz  -  1  kHz, 

1  Hz  -  10  kHz  and  1  Hz  -  25  kHz.  Since  the  measured  spectrum  is  only  an 
estimation  of  the  true  spectrum,  averaging  allows  us  to  improve  statistical 
accuracy  and  thus  to  reduce  the  variance  associated  with  the  measured 
spectrum.  So,  in  each  spectral  range  the  power  spectrum  obtained  was  the 
average  of  N  spectra  with  N  =  32.  The  results  reported  in  this  study  were 
collected  manually  and,  of  the  256  points  available  in  each  spectral  range, 
only  a  few  were  collected.  The  collected  points  were  chosen  so  that  they  would 
be  approximately  evenly  spaced  along  the  frequency  axis  on  a  logarithmic 
scale. 

For  ISFETs,  spectra  were  measured  in  1  mM,  10  mM,  and  100  mM  aqueous 
solutions  of  NaCl,  KC1  and  CaC^  for  Na+,  K+,  and  Ca^+  ISFETs,  respectively. 
For  pH  ISFETs,  spectra  were  measured  in  pH  4,  7,  and  10.  For  these 
measurements,  commercial  buffers  (VWR  Scientific  Inc.,  San  Francisco, 
California  94119)  were  used.  Finally,  for  redox  ISFETs,  spectra  were  measured 
for  different  ratios  of  [Fe(CN)^~]/[Fe(CN)^  ]  (0.1,  1,  and  10);  [Fe(CN)^-] 
being  kept  constant  at  10  mM;  these  solutions  were  prepared  using  K^Fe(CN)^ 
and  K^Fe(CN)^  and  were  buffered  at  pH  7.2  using  0.1  M  H^PO^.  All  chemicals 
used  were  of  reagent  grade.  Solutions  were  prepared  with  deionized  water 
having  a  resistivity  of  12  to  14  MO  cm.  All  measurements  were  conducted  at 
room  temperature. 

RESULTS  AND  DISCUSSION 
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Without  any  device  connected  in  the  system,  the  experimental  arrangement 

accurately  measures  the  Nyquist  noise  in  the  10  Kft  load  resistor,  as  shown  in 

Figure  3.  The  amplifier  noise  was  observed  to  be  significantly  lower 

-25  2 

(approximately  5  •  10  A  /Hz  for  10  Hz  <  f  <  25  kHz)  than  the  Nyquist  noise 

in  the  resistor  (the  amplifier  noise  represented  here  is  the  output  voltage 

2  8  2 

noise  spectral  density  divided  by  R*  =  10  ft  for  the  purpose  of  comparison 
with  1GFET  noise). 

Although  resistors  (especially  carbon  composition  resistors)  can  exhibit 
1/f  noise  when  current  flows  through  them  [19],  the  10  k£2  load  resistor  used 
did  not  show  noise  levels  larger  than  the  thermal  noise  when  a  current  of  a 
100  yA  was  passed  through  itf which  means  that  for  the  frequency  range  examined 
and  for  the  particular  level  of  current  chosen,  the  1/f  noise  is  hidden  by 
thermal  noise.  The  observed  noise  of  all  IGFETs  and  ISFETs  examined  was 
several  orders  of  magnitude  larger  than  the  amplifier  and  the  load  resistor 
noise  level. 

Figure  3  also  shows  the  current  noise  for  two  geometrically  identical 
IGFETs  processed  in  two  different  laboratories.  The  substantially  different 
noise  spectra  illustrate  the  dependence  of  the  "integrated  circuit"  noise  on 
the  processing  parameters  which  were  different  for  the  two  laboratories.  This 
variability  in  the  "integrated  circuit"  noise  was  present  to  a  much  smaller 
extent  In  devices  processed  in  the  same  laboratory  and  was  undetectable  for 
devices  from  the  same  wafer.  Spectra  (b)  and  (c)  in  Figure  3  are  most  likely  a 
superposition  of  1/f  noise  and  generation-recombination  (g-r)  noise  and 
similar  spectra  have  been  reported  in  the  literature  [20];  g-r  noise  may  be 
due  to  generation-recombination  centers  in  the  space  charge  region  existing 
between  the  conducting  channel  and  the  substrate;  on  the  other  hand,  1/f  noise 
may  be  due  to  interface  states  at  the  Si/SiOg  interface  and  to 


Craps  In  Che  SiOj  [12].  Spectrum  (c)  in  Figure  3  probably  indicaCes  a  higher 


concenCraCion  of  g-r  cenCers  in  Che  corresponding  IGFETs  Chan  in  chose 

exhibicing  spectrum  (b);  in  addition,  Che  noise  levels  reported  here  agree 

with  Chose  reported  for  metal  oxide  semiconductor  field  effect  transiscors 

(MOSFETs)  by  other  workers  [20,  21]. 

If  device  or  (Figure  1)  has  a  gold  gate,  it  can  be  operated  as  a 

regular  IGFET  by  directly  biasing  the  gate  or  it  can  be  operated  as  a  redox 

ISFET  by  biasing  the  gate  through  the  solution.  We  observed  no  difference  in 

the  noise  levels  between  these  two  modes  of  operation;  furthermore,  no 

4-  3- 

difference  was  observed  when  the  concentration  of  Fe(CN),  /Fe(CN),  redox 

o  o 

couple  was  changed  from  0.01  mM  to  10  mM  at  constant  ratio  of  1:1.  The 

standard  rate  constant  k°  for  the  charge  transfer  of  our  particular  couple 

—2  —1 

with  gold  is  on  the  order  of  10  cms  [22].  So,  this  change  in  concentration 

corresponds  to  a  change  in  exchange  current  density  j  at  the  gold/solution 

-5  2  -2  2 

interface  of  approximately  10  A/cm  to  10  A/cm  (Equation  9)  which  in  turn 

corresponds  to  a  variation  in  low  frequency  current  noise  spectral  density 
-20  2  -23  2 

of  about  3  •  10  A  /Hz  to  3  •  10  A  / Hz,  respectively  (Equations  10,  11, 

-4  2. 

and  15  and  having  a  gold  surface  area  of  2.3  *10  cm  ). 


j 


o 


i 

—2.  = 
A 


nFk° 


,  (l-a) 
0 


(9) 


where:  j  =  exchange  current  density 
iQ  =  exchange  current 
A  =  electrode  surface  area 
n  a  number  of  electrons  transferred 


F  =  Faraday  constant 


1 


k°  =  standard  rate  constant 
a  =  transfer  coefficient 

CQ  =  bulk  concentration  of  oxidized  species 
Cg  =  bulk  concentration  of  reduced  species 

The  relationship  between  exchange  current  and  noise  spectral  density 

will  be  discussed  below.  Considering  the  results  shown  in  Figure  4,  we  can 

immediately  see  that  these  values  are  below  the  "integrated  circuit"  noise. 

This  means  that  with  our  devices,  information  about  the  gold/solution 

interface  cannot  be  obtained  from  the  measurement  of  equilibrium  noise.  This, 

however,  is  a  technological  problem  and  commercially  available  field  effect 

-25  2 

transistors  can  reach  noise  levels  as  low  as  10  A  /Hz  (for  f  =  1  Hz)  which 
means  that  by  improved  processing  leading  to  a  decrease  in  "integrated 
circuit"  noise  we  may  indeed  be  able  to  get  information  about,  the 
gold/solution  interface  from  the  measurement  of  equilibrium  noise. 

An  interesting  difference  emerged  for  the  noise  spectrum  of  aluminum 
gate  IGFET  (Qj  or  Q^)  and  gold  gate  IGFET  (gate  or  covered  with  gold). 
The  latter  shows  higher  noise  levels  at  lower  frequencies  (Figure  4a)  than  the 
corresponding  aluminum  gate  IGFET  (Figure  4b)  on  the  same  chip.  This  observed 
additional  noise  may  be  explained  by  additional  interface  states  induced  by 
radiation  damage  during  sputtering  of  the  gold  gate  [23]  which  would  raise  the 
1/f  noise  level. 

The  interface  between  silicon  nitride  (Si^N^)  and  electrolyte  contains 
charge  which  is  pH  dependent  [ 1 ] .  Although  the  exact  nature  of  this  charge  is 
not  known,  ISFETs  with  Si^N^  are  used  as  pH  devices  which  are  free  from  common 
chemical  interference.  Their  response  to  steps  in  applied  gate  voltage  is  on 
the  order  of  a  microsecond  [24],  We  found  that  the  noise  spectrum  of  this 
device  is  undistinguishable  from  the  noise  of  the  corresponding  IGFET  with 
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aluminum  gate.  This  fact  is  in  agreement  with  our  interpretation  of  the  higher 
noise  levels  of  the  gold  gate  IGFETs  (Q^  or  Q^).  No  significant  difference  in 
the  noise  spectrum  of  pH  ISFETs  immersed  in  buffer  solutions  of  pH  4,  7,  and 
10  was  found.  As  in  the  case  of  redox  ISFE'rs,  it  can  be  concluded  that  the 
equilibrium  electrochemical  noise  is  obscured  by  the  "integrated-circuit" 
noise . 

A  different  situation  has  been  encountered  with  ISFETs  using  polymeric 
ion  selective  membranes.  Noise  spectra  of  Na+  and  K+  ISFETs  are  shown  in 
Figure  5,  together  with  the  spectrum  of  pH  ISFETs.  The  polymeric  membranes  are 
hydrated  which  means  that  the  silicon  nitride  surface  under  the  membrane  is 
also  hydrated.  From  this  point  of  view  it  is  appropriate  to  subtract  the  noise 
spectrum  of  the  pH  ISFET  from  the  noise  spectra  of  polymeric  membrane  ISFETs. 
Spectra  for  Na+  and  K+  ISFETs  in  Figure  4  were  obtained  after  subtracting  the 
pH  ISFET  spectrum.  For  frequencies  higher  than  100  Hz,  K+  ISFETs  did  not  show 
noise  levels  significantly  higher  than  those  of  pH  ISFETs,  which  is  the  reason 
why  only  points  below  100  Hz  are  shown  in  this  case.  Since  devices  are 
operated  in  saturation,  current  in  the  10  kfi  load  resistor  is  effectively  the 
short  circuit  drain  current.  Although  the  noise  spectral  density  is  elevated 
by  the  presence  of  polymeric  ion  selective  membranes,  no  significant  variation 
was  observed  for  different  concentrations  of  the  primary  ion  (1,  10,  100  mM). 
The  noise  spectrum  of  Ca^+  ISFETs  was  only  slightly  higher  that  that  of  the  pH 
ISFET  for  frequencies  larger  than  100  Hz  and  was  not  included  in  Figure  5. 

Impedance  measurements  applied  to  the  investigation  of  ion  selective 
electrodes  [25,  26]  show  that  such  electrochemical  systems  can  be  modeled  with 
equivalent  electrical  circuits  composed  of  parallel  RC  combinations  in  series, 
each  combination  corresponding  to  a  specific  characteristic  of  the  system. 

This  approach  is  very  Important  for  a  better  understanding  of  such  systems  and 
is  useful  in  the  identification  of  the  parameters  limiting  the  performance  of 
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these  systems.  In  thermodynamic  equilibrium,  every  impedance  Z  exhibits 
thermal  noise,  the  voltage  spectral  density  S^C f )  of  which  is  given  by  the 
Nyquist  equation: 


Sy(f)  =  4  kT  Re(Z) 


(10) 


where  k  =  Boltzmann  constant 

T  =  absolute  temperature 

Re(Z)  =  real  part  of  impedance  Z 

This  is,  of  course,  white  noise. 

Using  computer  assisted  modeling,  we  attempted  to  determine  an 
equivalent  electrical  circuit  for  Na+  and  K+  ISFETs  so  that  the  theoretical 
noise  spectra  would  fit  the  experimental  points  as  shown  in  Figure  5.  For  this 

task,  we  used  a  general  purpose  circuit  simulation  program  (SPICE)  developed 

by  the  Computer  Sciences  Department  at  the  University  of  California,  Berkeley 
and  adapted  for  use  on  a  HP  3000  computer  by  the  Hewlett-Packard  Design  Aids 
Group  [27].  Figure  6a  shows  the  electrical  equivalent  circuit  used  for 
modeling.  This  circuit  is  comprised  of  two  parts.  The  first  part  models  the 
electrochemical  system  placed  over  the  integrated  circuit  and  is  made  up  of  a 


series  of  parallel  R^C^  combinations;  at  this  point  we  do  not  assign 
these  Rj  and  elements  to  any  physical  part  of  the  electrochemical  system; 
the  resistances  R^  are  generating  thermal  noise  at  thermal  equilibrium; 
in  our  model  we  represent  them  as  noiseless  elements  and,  in  series  with  them, 


we  place  a  noise  source  which  has  a  spectral  density  given  by  Equation  10 
and  which  is  equivalent  to  the  thermal  noise  generated  across  Rj.  The  second 
part  of  the  circuit  models  the  integrated  circuit  and  is  the  simplest 
equivalent  circuit  for  an  1GFET  in  saturation  [28].  In  this  model  C__ 
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represents  the  gate  to  source  capacitance,  g  is  the  transconductance  of  the 
■■  m 

device,  V^g  is  the  voltage  developed  across  C^g,  1Q  is  the  drain  current  and 

is  equal  to  g  V  and  finally  R.  is  the  load  resistance  (equivalent  to  R.  in 
ro  ub  ti  4 

-w  —U  —1 

Figure  2).  For  our  devices,  we  had  C-c  =  1.5  pF  and  g  =  4  •  10  ft  . 

Ob  m 

It  can  be  safely  assumed  that  the  noise  sources  V  due  to  the 

1/2 

resistances  R,  are  uncorrelated.  V.  is  such  that  V.  =  (S  .(f))  and  because 
j  J  J  vj 

noise  adds  by  adding  spectral  densities,  we  can  write  for  the  spectral 
density  S^(f)  of  the  drain  currents 

Si<£>  -  8„2  5  Vj2  <u> 

This  relationship  is  exact  only  at  zero  frequency  but  it  applies  quite  well 

for  low  frequencies  at  which  filtering  due  to  capacitances  C,  and  C  is 

J  Ob 

negligible.  Figure  6b  schematically  shows  how  the  experimental  spectrum  S^f) 
is  taken  as  being  the  sum  of  several  contributions  S^(f)  due  to  the  thermal 
noise  in  R ^ .  The  mathematical  expression  for  the  frequency  dependence  of 
S ^ j ( f )  is  very  complicated  and  will  not  be  given  here.  The  values  for  R^  were 
obtained  by  using  the  expression: 

S. .(0) 

R.  =  - U—T  (12) 

J  4  kT  g  * 

Ttl 

which  is  a  combination  of  Equations  10  and  11.  On  the  other  hand  initial 
values  were  obtained  using  the  approximate  expression: 

sn<0) 

S  <£>  ; - 5-5±J - 2  (13) 

1  1  ♦  “  <cj  *  CGS> 

This  expression  allows  us  to  obtain  approximate  values  for  Cj  from  the 
experimental  spectrum.  These  values  are  then  varied  in  order  to  improve  the 


fit  of  the  model  with  the  experimental  data.  Therefore,  using  Equation  13  and 
obtaining  the  half-power  frequency  f ^  from  the  experimental  spectrum  (Figure 


6b)  we  can  write: 


C j  “  2  TT  f  ~  CGS 


(14) 


The  modeling  for  a  given  device  is  done  with  a  minimal  number  of  RC 
combinations  required  to  obtain  a  satisfactory  fit.  The  values  which  were  used 
to  fit  the  data  presented  in  Figure  5  are  shown  in  Table  I.  The  correlation 
coefficient  r  between  the  model  and  *he  experimental  data  was  calculated  and 
we  obtained  r  =  0.9782  for  the  Na+  ISFET  and  r  =  0.9980  for  the  K+  ISFET.  For 
the  K+  ISFET,  the  correlation  coefficient  was  calculated  for  5  <  f  <  100  Hz  , 
since  for  frequencies  higher  than  100  Hz  the  noise  spectral  density  for  this 
device  does  not  differ  from  the  one  for  pH  ISFETs.  As  shown  in  Table  I,  we 
have  modeled  the  Na+  ISFET  with  three  RC  combinations;  only  one  RC  combination 
was  used  for  the  small  frequency  range  in  which  the  K+  ISFET  can  be  modeled. 
Since  we  observe  no  changes  in  noise  spectral  density  when  changing  the  main 
ion  concentration  in  solution,  it  must  be  that  the  dominant  noise  is  due  to  a 
resistance  which  does  not  change  its  value  under  these  conditions.  Such  a 
resistance  may  well  be  the  bulk  resistance  of  the  membrane  and  so,  the  R^C^ 
combination  may  be  assigned  to  the  bulk  properties  of  the  membrane;  this 
applies  for  both  K+  and  Na+  ISFETs.  The  R^C^  combination  observed  for  the  Na+ 
ISFETs  is  most  likely  associated  with  the  membrane/silicon  nitride  interface 
since  the  time  constant  Rj(Cj  +  C^g)  "  ^  is  on*y  an  or^er  magnitude 
larger  than  the  electrical  response  time  of  the  pH  ISFET  which  has  been 
measured  to  be  on  the  order  of  a  microsecond  [24].  Finally,  RjCj  he  due  to 
the  solution/membrane  interface  and  thus,  Rj  would  be  the  charge  transfer 
resistance  R  and  would  be  the  double  layer  capacitance  of  this  particular 
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TT 


interface.  The  charge  transfer  resistance  R  is  inversely  proportional  to  the 

ct 

exchange  current  iQ,  as  shown  by  Equation  15: 


R 


ct 


RT 

nFi 


(15) 


where  R  =  Molar  gas  constant 

T  =  absolute  temperature 

n  =  number  of  elementary  charges  transferred  across  the  interface  per 
elementary  event 
F  =  Faraday  constant 

iQ  =  exchange  current  across  the  interface 
Rct  should  thus  change  with  varation  of  the  exchange  current  i^  which  would 
lead  to  a  change  in  the  corresponding  noise  spectral  density;  according  to 
Equation  9,  iQ  should  change  with  the  ionic  activity  of  the  exchanging  ion. 
However,  we  have  mentioned  earlier  that  no  change  in  noise  spectral  density 
was  observed  when  varying  the  solution  ionic  activity.  It  is  quite  conceivable 
that  Rct  indeed  changes  with  varying  ionic  activity  but  that  this  change  is 
obscured  by  the  noise  of  larger  magnitude  due  to  the  bulk  resistance  of  the 
membrane.  For  the  Na+  ISFET  the  value  of  Rj  used  in  Equation  15  leads  to  an 
exchange  current  across  the  interface  of  approximately  2  10  10  A.  The 

polymeric  membranes  which  were  used  in  this  study  had  poorly  defined  membrane 
surface  area  which  makes  the  interpretation  of  our  results  in  terms  of 

2 

exchange  current  density  difficult;  however,  if  we  consider  an  area  of  1  mm 

which  is  a  reasonable  estimate,  it  leads  to  an  exchange  current  density  of  2 
—8  2 

10  A/cm  ,  which  is  a  rather  small  value  for  a  nonpolarized  interface, 
considering  the  solution  ionic  activities  used  in  our  experiments.  A  study  on 
Na+  ISFETs  [18]  reported  an  electrical  response  time  (response  to  a  voltage 
step  applied  to  the  reference  electrode)  of  approximately  15  ms;  T  is  the  time 
required  to  reach  63%  of  the  full  response.  If  we  apply  a  voltage  step  to 
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our  model  for  the  Na+  ISFET,  we  obtain  a  value  of  about  0.2  ms  for  f  This 
apparent  discrepancy  can  be  explained  by  considering  the  simple  circuit  shown 
in  Figure  7a.  R  and  C  represent  the  impedance  of  a  hypothetical  membrane  and 
Cgg  represents  the  input  capacitance  of  our  integrated  circuit,  when  operated 
in  the  saturation  mode.  If  V^(t)  is  a  step  function  of  amplitude  VQ,  the  time 
course  of  the  voltage  VGg(t)  across  CGg  is  given  by  Equation  16: 


V0SM  -  vo  (1  ' 


e-t/E(C  .  CGS>) 


C  -f  c. 


For  t  =  0,  VQg(0)  is  given  by  Equation  17, 


VGS(0)  =  V0  C  +  C, 


This  initial  fast  step  is  always  present  in  our  situation;  however,  it  is  only 
apparent  for  C  >  0.1  CGg.  After  the  initial  fast  step  is  reached,  the  response 
reaches  exponentially  with  a  time  constant  T*  equal  to  R(C  +  CGg).  If  R  and 
C  are  due  only  to  the  geometry  of  the  membrane,  even  though  R  and  C  are 
geometry  dependent,  the  product  RC  is  indeed  geometry  independent.  However, 
because  of  the  low  CGg  values  encountered  in  ISFETs,  the  condition  C  0.1  CGg 
is  not  always  satisfied  and  a  situation  similar  to  the  one  illustrated  in 
Figure  7b  can  arise.  The  simulation  was  done  for  a  constant  RC  product  with 
the  ratio  a  =  C/C„_  as  a  parameter.  Figure  7b  shows  the  effect  of  membrane 
geometry  on  the  electrical  response  time  of  the  whole  system.  Returning  to  the 
Na+  device,  we  can  see  that  the  initial  fast  step  is  determined  by  the 
capacitive  divider  formed  by  CGg  and  the  series  combination  of  C^,  ,  and 

(Figure  6a)  which  according  to  the  values  shown  in  Table  1  would  correspond  to 
an  initial  step  amounting  to  approximately  U2X  of  the  full  response  and  would 
give  a  response  shape  similar  to  the  one  for  a  =  1  in  Figure  7b.  Because  of 


differences  in  the  membrane  geometry,  it  may  well  be  that  the  results 
presented  in  Reference  13  have  a  response  shape  intermediate  between  those  for 
a  =  0.2  and  a  =  0.01  in  Figure  7b, and  thus  would  show  a  T  significantly  larger 
than  the  value  presented  here.  In  the  case  of  the  K+  ISFET,  we  cannot  talk 
about  response  time  values  since  ^2^2  an<*  R3^3  cou^  not  he  identified  even 
though  they  are  physically  present  in  the  device,  and  as  shown  in  the  case  of 
the  Na+  ISFET,  will  affect  the  shape  of  the  response  curve. 

It  has  been  pointed  out  [5]  that  equilibrium  noise  measurements  can  give 
information  about  the  frequency  dispersion  of  the  resistive  and  reactive  parts 
of  the  equivalent  electrical  circuit  for  a  system  as  well  as  information  about 
the  exchange  current  density  at  the  membrane/ solution  interface.  For  metal 
electrodes,  this  value  can  be  obtained  from  Tafel  plots,  however,  for  ion 
selective  membranes  in  which  the  bulk  resistance  is  comparable  to  the  charge 
transfer  resistance,  Tafel  plots  cannot  be  obtained  and  the  measurement  of 
equilibrium  noise  could  be  the  only  way  to  obtain  this  information.  The  use  of 
ISFETs  is  justified  because  of  the  wide  bandwidth  characteristics  of  these 
devices  and  the  small  surface  area  of  their  ion  selective  membrane.  A  better 
definition  of  membrane  geometrical  characteristics  can  be  achieved  by  using 
suspended  mesh  ISFETs  which  were  recently  developed  [29], 

The  present  study  has  shown  two  important  limitations  of  the  technique. 
First  of  all,  the  "integrated  circuit"  noise  depends  strongly  on  processing. 

It  is  not  possible  to  compare  directly  the  noise  spectra  of  devices  obtained 
from  different  laboratories.  Secondly,  the  "integrated  circuit"  noise  imposes 
a  background  level  which  must  be  exceeded  by  the  noise  of  the  electrochemical 
system  placed  on  the  integrated  circuit  in  order  to  extract  any  information 
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from  it.  For  example,  the  noise  spectra  for  Ca^+  and  to  some  extent  K+  ISFETs 
were  too  close  to  the  noise  spectrum  of  the  corresponding  pH  ISFET  in  order  to 
obtain  kinetic  information  in  the  frequency  range  investigated.  Since  most  ion 
selective  membranes  and  other  nonpolarized  electrodes  have  exchange  current 
densities  higher  than  the  values  found  in  this  study  for  the  Na+  ISFET,  they 
probably  could  not  be  investigated  with  this  method  if  the  "integrated 
circuit"  noise  is  not  reduced  by  improved  processing.  It  is,  however,  possible 
that  noise  measurements  made  on  these  systems  while  they  are  driven  out  of 
equilibrium  (non-zero  net  current)  will  provide  additional  information  about 
the  processes  involved  at  the  solution/membrane  interface. 
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FIGURE  CAPTIONS 


Figure  1.  Photograph  of  integrated  circuit  showing  the  two  ISFETs  (Q  ,  Q^) 
and  the  two  IGFETs  (Qj,  Q^). 


Figure  2.  Experimental  arrangement  for  ISFET  and  IGFET  noise  measurements 

(shown  with  IGFET).  Values  are:  Rx  =  25  kfl  ±  1%;  R2  =  0  -  50  kft 
(potentiometer);  R^  =  100  kfi  -  1*4;  R^  =  10  kf2  ±  1*4;  C  =  1  yF;  = 
9  V;  V2  =  6  V. 


Figure  3.  Noise  Spectra. 


a)  Q]  :  amplifier  +  resistor  noise. 

b)  0  :  IGFET  processed  in  the  HEDCO  Laboratory  of  the  University 

of  Utah  (average  of  five  measurements). 

c)  •  :  IGFET  processed  in  the  microcircuit  laboratory  of  the 

University  of  Utah  Research  Institute  (Average  of  five 
measurements) . 

The  dotted  line  indicates  the  theoretical  level  for  thermal  noise 
in  the  10  resistor.  (This  level  is  1.66  •  10  ^  A^/Hz  at  300°K> 


Figure  4.  Noise  spectra. 


a)  •  :  gold  gate  IGFET  (Q^  or  Q^)  (average  of  five  measurements) 
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b)  0  :  aluminum  gate  IGFET  (Qj  or  Q^)  on  same  substrate  (average 

of  five  measurements) 


Figure  5.  Noise  Spectra. 

a)  0  :  pH  ISFET  (average  of  five  measurements) 

b)  •  s  K+  ISFET  (average  of  three  measurements)  pH  noise  has  been 
subtracted  (dotted  lines  corresponds  to  model) 

<=>□  :  Na+  ISFET  (average  of  three  measurements)  pH  noise  has 
been  subtracted  (solid  line  corresponds  to  model) 

Figure  6.  a)  Electrical  equivalent  circuit  used  for  modeling. 

b)  schematic  explanation  of  how  the  experimental  spectrum  is  fit 
with  a  sum  of  individual  spectra,  each  corresponding  to  a  single 
noise  source  V^.  (A  log-log  representation  is  shown). 

Figure  7.  a)  Circuit  showing  a  hypothetical  membrane  having  an  impedance 
equivalent  to  R  and  C  in  parallel;  CQg  represents  the  input 
capacitance  of  the  integrated  circuit. 

b)  Computer  simulated  response  for  circuit  shown  in 

part  (a)  for  being  a  voltage  step.  Simulation  is  done  for  RC  = 
constant  =  1  ms  and  for  various  ratios  o(  a  <  C/C__;  (C__  * 

VlS  V»9 

constant  *  1  pF).  T  and  T*  are  shown  for  a  *  0.2.  T  is  measured  for 
the  response  from  the  baseline  and  T*  is  measured  for  the  response 
from  the  initial  fast  step. 
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TABLE  I 


Computer  Modeling  Results  for  Na+  and  K+  ISFETs 


1SFET 

m 

C1 

[pF] 

m 

u 

l^j 

C3 

tpF] 

Na+ 

1130.0 

O 

113.0 

D 

3.8 

5.5 

K+ 

189.0 

26.6 

— 

— 

— 

1 

1 

=  membrane  bulk  properties 
R2C2  £  membrane/solution  interface 
R^C^  =  membrane /Si^N^  interface 
(See  text) 
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